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C^ 
'~H ^ Quantum plasmas is a rapidly expanding field of research, with applications 

5-^1 ranging from nanoelectronics, nanoscale devices and ultracold plasmas, to in- 

ly^ ' ertial confinement fusion and astrophysics. Here we give a short systematic 

O ' overview of quantum plasmas. In particular, we analyze the collective effects 

' J-»^ , due to spin using fluid models. The introduction of an intrinsic magnetization 

^. ■ due to the plasma electron (or positron) spin properties in the magnetohydro- 

(^ ' dynamic limit is discussed. Finally, a discussion of the theory and examples of 

f^ ' applications is given. 

1. Introduction 

_^ ■ The field of quantum plasmas is a rapidly growing field of research. From 

fT^ , the non-relativistic domain, with its basic description in terms of the 

C^ ' Schrodinger equation, to the strongly relativistic regime, with its natu- 

ral connection to quantum field theory, quantum plasma physics provides 
'nI I promises of highly interesting and important application, fundamental con- 

^"^ ■ nections between different areas of science, as well as difficult challenges 

f^ , from a computational perspective. The necessity to thoroughly understand 

such plasmas motivates a reductive principle of research, for which we suc- 
cessively build more complex models based on previous results. The simplest 
. , lower order effect due to relativistic quantum mechanics is the introduction 

d ' of spin, and as such thus provides a first step towards a partial description 

of relativistic quantum plasmas. 

Already in the 1960's, Pines studied the excitation spectrum of quan- 
tum plasmas, ^'^ for which we have a high density and a low temperature as 
compared to normal plasmas. In such systems, the finite width of the elec- 
tron wave function makes quantum tunnelling effects crucial, leading to an 
altered dispersion relation. Since the pioneering work by Pines, a number 
or theoretical studies of quantum statistical properties of plasmas has been 
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done (see, e.g., Ref. 3 and references therein). For example, Bezzerides 
& DuBois presented a kinetic theory for the quantum electrodynamical 
properties of nonthermal plasmas,^ while Hakim & Heyvaerts presented 
a covariant Wigner function approach for relativistic quantum plasmas.^ 
Recently there has been an increased interest in the properties of quan- 
tum plasmas. ®^^^ The studies has been motivated by the development in 
nanostructured materials^^ and quantum wells, ^^ the discovery of ultra- 
cold plasmas^^ (see Ref. 29 for an experimental demonstration of quantum 
plasma oscillations in Rydberg systems), astrophysical applications,'^*' or a 
general theoretical interest. Moreover, it has recently been experimentally 
shown that quantum dispersive effects are important in incrtial confinement 
plasmas. '^^ The list of quantum mechanical effects that can be included in a 
fluid picture includes the dispersive particle properties accounted for by the 
Bohm potential,®"^* the zero temperature Fermi pressure, ^"^^ spin proper- 
ties^^"^^ as well as certain quantum electrodynamical effects. ^^"^''' Within 
such descriptions, ^"^"'^^'^"'^'^"'^^ quantum and classical collective effects can 
be described within a unified picture. 

2. The microscopic equations: Schrodinger and Pauli 
dynamics 

2.1. The Schrodinger description 

The basic equation of nonrelativistic quantum mechanics is the Schrodinger 
equation. The dynamics of an electron, represented by its wave function 2p, 
in an external electromagnetic potential (j) is governed by 

where h is Planck's constant, me is the electron mass, and e is the mag- 
nitude of the electron charge. This complex equation may be written as 
two real equations, writing tp = ^fn exp iSjh, where n is the amplitude 
and 5* the phase of the wave function, respectively.^® Such a decomposition 
was presented by de Broglie and Bohm in order to understand the dynam- 
ics of the electron wave packet in terms of classical variables. Using this 
decomposition in Eq. (1), we obtain 

^ + V.(nv)^0, (2) 

and 

me— = eV0+^^V^^^ , 3 

at 2me \ Jn / 



April 2, 2008 11:40 WSPC - Proceedings Trim Size: 9in x 6in Trieste07-spin-MHD 



where the velocity is defined by v = VS/nie. The last term of Eq. (3) is 
the gradient of the Bohm-de Broglie potential, and is due to the effect of 
wave function spreading, giving rise to a dispersive-like term. We also note 
the striking resemblance of Eqs. (2) and (3) to the classical fluid equations. 

2.2. The Pauli description 

In relativistic quantum mechanics, the spin of the electron (and positron) 
is rigorously introduced through the Dirac Hamiltonian 

H ^ ca ■ {p + eA) - ecj) + (3meC^ , (4) 

where a = (ai, a2, as), e is the magnitude of the electron charge, c is the 
speed of light, A is the vector potential, (j) is the electrostatic potential, and 
the relevant matrices are given by 

Here I is the unit 2x2 matrix and a = (cri, (72, era), where we have the Pauli 
spin matrices 

"^^(iJ)'"'=(? o')' ""^"^=((l-i)- (^^ 

From the Hamiltonian (4), a nonrclativistic counterpart may be obtained, 
taking the form 

H^J-ip + eAf + ^B-a-ec^. (7) 

Thus, the electron possesses a magnetic moment m — — /iB(V'l'''l'0)/('0IV'): 
where /is = eh/2me is the Bohr magneton, giving a contribution — B • m 
to the energy. The latter shows the paramagnetic property of the electron, 
where the spin vector is anti-parallel to the magnetic field in order to mini- 
mize the energy of the magnetized system. According to (7) and the relation 
dF/dt = dF/dt + {l/iH)[F, H], where F is some operator and [, ] is the Pois- 
son bracket, we have the following evolution equations for the position and 
momentum in the Hcisenberg picture^®"'^* 

me-^ = -e(E + vxB)--MsV(B-s), (9) 



April 2, 2008 11:40 WSPC - Proceedings Trim Size: 9in x 6in Trieste07-spin-MHD 



while the spin evolution is given by 



ds 2 ^ 

— = -mbB X s, (10) 

where the spin operator is given by 

s^-a. (11) 

The above equations thus gives the quantum operator equivalents of the 
equations of motion for a classical particle, including the evolution of the 
spin in a magnetic field. 

The non-relativistic evolution of spin i particles, as described by the 
two-component spinor ^(a), is given by the Pauli equation (see, e.g., 36) 

2 " 






2me \ h 



ZG 

V + —A I - /isB -a + ecf) 



V'-O, (12) 



where A is the vector potential, jib = eh/2mf, is the Bohr magneton, and 
a = ((Ti, (72, (Ta) is the Pauli spin vector. 

Now, in the same way as in the Schrodingcr case, we may decompose 
the electron wave function tp into its amplitude and phase. However, as 
the electron has spin, the wave function is now represented by a 2-spinor 
instead of a c-number. Thus, we may use -0 — ^/n exp(iS'/ft)(/3, where (p, 
normalized such that (p^ip = 1, now gives the spin part of the wave function. 
Multiplying the Pauli equation (12) by tp\ inserting the above wave func- 
tion decomposition and taking the gradient of the resulting phase evolution 
equation, we obtain the conservation equations 

^+V.(nv)=0 (13) 

and 

dv /^ x.^ ^^ - /V^ 

TOe — = — e(E + V X B) 
dt 



ft2 


YV 


2m, 






1 



-^(V®B)-s-^V-(nE) (14) 

h rrien 

respectively. The spin contribution to Eq. (14) is consistent with the results 
of Ref. 39. Here the velocity is defined by 

v^ — (VS-ihtp''Vip) + —, (15) 

nie nieC 



the spin density vector is 

h 



^ip^'aif, (16) 
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which is normalized according to 

|s| = h/2, (17) 

and we have defined the symmetric gradient spin tensor 

Z = (Vs,)®(Vs''). (18) 

Moreover, contracting Eq. (12) by i/i^cr, we obtain the spin evolution equa- 
tion 

§ = (^B-^[a„(na"s)]jxs. (19) 

(it [^ ft m^n J 

We note that the last equation allows for the introduction of an effective 
magnetic field Bcff = (2/iB/fi,)B — {nien)^^ [da{nd°'s)\. However, this will 
not pursued further here (for a discussion, see Ref. 36). 

Comparing the effects due to spin from the Pauli dynamics with the 
Schrodinger theory, we see a significant increase in the complexity of the 
fiuid like equations due the presence of spin. The fact that the spin couples 
linearly to the magnetic field makes the dynamical aspects of such Pauli sys- 
tems very rich. Moreover, when going over to the collective regime, the back 
reaction through Maxwell's equation can yield interesting new properties 
of such spin plasmas. In fact, the introduction of an intrinsic magnetization 
can give rise to linear instability regimes, much like the Jeans instability 
(see Sec. 5.2.). 

3. Collective plasma dynamics 

As pointed out in the previous section, the route from single wavefunction 
dynamics to collective effects introduces a new complexity into the system. 
At the classical level, the ordinary pressure is such an effect. In the quantum 
case, a similar term, based on the thermal distribution of spins, will be 
introduced. 

3.1. Multistream model 

The multistream model of classical plasmas was successfully introduced 
by Dawson. *° Here we will focus on the electrostatic interaction between 
a multistream quantum plasma described within the Schrodiner model, a 
system first investigated in Ref. 7 (where also the stationary regime was 
probed). Thus, we have the governing equations (2) and (3) but for N beams 
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of electrons on a stationary ion background, i.e., Using this decomposition 
in Eq. (1), we obtain 

^+V-(n„v„)=0, (20) 



and 



dt 



m^^=eW<P+^w(^^^=^], (21) 



dt 2me 

now coupled through the self-consistent electrostatic potential governed by 

N 

V'<^=- V(n„-no). (22) 

eo 1 

a— 1 

Here, d/dt = 9^ + v^ • V and no is the density of the stationary ion back- 
ground. 

In the one-stream case {a — 1), we have the equilibrium solution v = vq 
(a constant drift relative the stationary ion background) and the constant 
electron density n = no (such that = 0). Perturbing this system a Fourier 
decomposing the perturbations, such that n — ng + (5nexp[z(k • x — cot)]), 
V = Vo -I- (5vexp[i(k • x — cot)], and cj) = S<j)exp[i{'k ■ x — cot)], we obtain-^'"^ 

4to2 ' 



{u-k- vo)2 = ^l + —^, (23) 



where the last term is the Bohm-de Broglie correction to the dispersion 
relation. Here we have the electron plasma frequency ujp = (e^no/eome)^'^. 
Similarly to the one-stream case, we obtain the dispersion relation^''*^ 



(W - vox • k)2 - ft2fc4/4. 



TO^ 



+ ^ (24) 

(w - V02 • k)2 - K^k^/Aml ' ^ ' 

for two propagating electron beams (with velocities vqi and V02) with back- 
ground densities noi and no2. The quantum effect has a subtle influence on 
the stability of the perturbed plasma. For the case noi = rio2 = no/2 and 
Vqi = — V02 = vq, we have the instability condition 

^(l-^)<^^<^^ (25) 

in terms of the normalized wavenumber K = kvo/ujp and the quantum 
parameter H = hajp/nieVQ (see Fig. l).^'**^ We see that when H ~ 0, we 
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Stable 
region 




Vok/Ojp 

Fig. 1. The regions of stability and instability in the case of the quantum two-stream 
interaction.^'*^ 

have unstable perturbations for < X < 1, but when H y^ a, considerably 
more complex instability region develops. 

A model for treating partial coherence in such systems, based on the 
Wigner transform technique, ''^"''^ can also be developed'*^ (see also Ref. 
46). Moreover, using the equations (13) and (14), a similar framework may 
be set up for electron streams with spin properties. 



3.2. Fluid model 

3.2.1. Plasmas based on the Schrodinger model 

Suppose that we have N electron wavef unctions, and that the total system 
wave function can be described by the factorization '0(xi,X2, . . .x^v) = 
V'iV'2 ■ • ■ "fpN- For each wave function ipa, we have a corresponding proba- 
bility Va- From this, we first define "00. — naexp{iSa/fi) and follow the 
steps leading to Eqs. (2) and (3). We now have N such equations the wave 
functions {ipa}- Defining^ 



JV 

E 



/ a^a 



(26) 
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and 



N „ 

J2 ^"^"^" , (27) 

a— 1 



we can define the deviation from tlie mean flow according to 

Wq =: v„ - V. (28) 

Taking the average, as defined by (27), of Eqs. (2) and (3) and using the 
above variables, we obtain the quantum fluid equation 

^+V-(nv)=0 (29) 



and 



men ( — + V • V j V = enVcj) ~Vp+ ^V / ( — ^ ) ) , (30) 



where we have assumed that the average produces an isotropic pressure 
p = TOen(|wQ,p) We note that the above equations still contain an explicit 
sum over the electron wave functions. For typical scale lengths larger than 
the Fermi wavelength Xp, we may approximate the last term by the Bohm~ 
de Broglie potential^ 



(31) 

Using a classical or quantum model for the pressure term, we finally have 
a quantum fiuid system of equations. For a self-consistent potential (f> we 
furthermore have 

V^(j)^ —{n-Ui). (32) 

eo 

3.2.2. Spin plasmas 

The collective dynamics of electrons with spin and some of the spin modi- 
fications of the classical dispersion relation was presented in Ref. 19. Here 
we will follow Refs. 19 and 20 for the derivation of the governing equations. 
Suppose that we have N wave functions for the electrons with magnetic 
moment /Xg = —pLsi a-nd that, as in the case of the Schrodinger descrip- 
tion, the total system wave function can be described by the factorization 
■)/) = ■0]^-(/'2 . . -ipN- Then the density is defined as in Eq. (26) and the aver- 
age fluid velocity defined by (27). However, we now have one further fluid 
variable, the spin vector, and accordingly we let S = (s^). From this we 
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can define the microscopic microscopic spin density <Sq = Sq, — S, such that 
(«S„) = 0. 

Taking the ensemble average of Eqs. (13) we obtain the continuity equa- 
tion (29), while we the the ensemble average applied to (14) yield 



— + V • V V = -en (E + V X B) - Vp + V — ^ 

ot J 2me \ vn 

(33) 
and the average of Eq. (19) gives 



I— )^^T' 



n ( — + V • V ) S = -4^B X S - V • K + fispin (34) 



respectively. Here the force density due to the electron spin is 

F,pi„ ^ -^(V ® B) . S - -i-V • [n{-L + T. ] 

-—V • [n{VSa) (E> (V5^) + n{VS^) (g> (V5„)] , (35) 
nie 

consistent with the results in Ref. 39, while the asymmetric thermal-spin 
coupling is 

K = n{wa(dSa) (36) 

and the nonlinear spin fluid correction is 

apin = — S X [daind'^S)] + — S X [da{n{d''Sa.))] 

rrie me 

+ ^/^x{a„Ka"(S + 5„)]}\ (37) 

nie \na I 

where Z = i^Sa) <E) (VS") is the nonlinear spin correction to the clas- 
sical momentum equation, Z = {{'^S(a)a) (^ (^^(a))) ^^ ^ pressure like 
spin term (which may be decomposed into trace- free part and trace), and 
[(V B) • S]" = {d''Bb)SK Here the indices a,6, . . . = 1, 2,3 denotes the 
Cartesian components of the corresponding tensor. We note that, apart 
from the additional spin density evolution equation (34), the momentum 
conservation equation (33) is considerably more complicated compared to 
the Schrodinger case represented by (30). Moreover, Eqs. (33) and (34) still 
contains the explicit sum over the N states, and has to be approximated 
using insights from quantum kinetic theory or some effective theory. 

The coupling between the quantum plasma species is mediated by the 
electromagnetic field. By definition, we let H = B//xo — M where M = 
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— 2n/iBS/fi, is the magnetization due to the spin sources. Ampere's law 
V X H = j + eo^iE takes the form 

1 fiV 

VxB = /io(j+VxM) + - — , (38) 

c' at 

where j is the free current contribution The system is closed by Faraday's 

law 

VxE^--. (39) 

4. The magnetohydrodynamic limit 

The concept of a magnetoplasma was first introduced in the pioneering work 
47 by Alfven, who showed the existence of waves in magnetized plasmas. 
Since then, magnetohydrodynamics (MHD) has found applications in a vast 
range of fields, from solar physics and astrophysical dynamos, to fusion 
plasmas and dusty laboratory plasmas. 

Magnetic fields, an essential component in the MHD description of plas- 
mas, also couples directly to the spin of the electron. Thus, the presence 
of spin alters the single electron dynamics, introducing a correction to the 
Lorentz force term. Indeed, from the experimental perspective, a certain 
interest has been directed towards the relation of spin properties to the 
classical theory of motion (see, e.g., Refs. 48-60). In particular, the effects 
of strong fields on single particles with spin has attracted experimental 
interest in the laser community. '''""^^ However, the main objective of these 
studies was single particle dynamics, relevant for dilute laboratory systems, 
whereas our focus will be on collective effects. 

We will now include if the ion species, which are assumed to be described 
by the classical equations and have charge Ze, we may derive a set of one- 
fluid equations.'^" The ion equations read 

^+V.(n,v,) = 0, (40) 

at 

and 

d 



Q^ , V, • V J V, = Zen, (E + v, x B) - -Vp,. (41) 

Next we define the total mass density p = {men + rmni), the centre-of- 
mass fiuid flow velocity V = [nienWf, +miniWi)/ p, and the current density 
j ~ — enve + Zeui-Vi. Using these denfinitions, we immediately obtain 

^+V.(pY)^0, (42) 
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from Eqs. (29) and (40). Assuming quasi-neutrality, i.e. n sa Zrii, the mo- 
mentum conservation equations (33) and (41) give 

| + v-v)v.jxB-v.n-v,.|gLv(^)+F.,„.(«, 

where fl is the tracefree pressure tensor in the centre-of-mass frame, and 
P is the scalar pressure in the centrc-of-mass frame. We also note that due 
to quasi-neutrality, we have Ug ~ Zp/rrii and v = V — mi']/ Zep, and we 
can thus express the quantum terms in terms of the total mass density p, 
the centre-of-mass fluid velocity V, and the current j. With this, the spin 
transport equation (34) reads 

In the momentum equation (43), neglecting the pressure and the Bohm- 
de Broglie potential for the sake of clarity, we have the force density j x 
B-|-Fspin. In general, for a magnetized medium with magnetization density 
M, Ampere's law gives the free current in a finite volume V according to 

j = —V X B - V X M, (45) 

where we have neglected the displacement current. The surface current is 
an important part of the total current when we are interested in the forces 
on a finite volume, as was demonstrated in Ref. 20 and will be shown below. 
It it worth noting that the expression of the force density in the momen- 
tum conservation equation can, to lowest order in the spin, be derived on 
general macroscopic grounds. Formally, the total force density on a volume 
element V is defined as F = \m\v^[){Y^^ ia/V), where f^ are the different 
forces acting on the volume element, and might include surface forces as 
well. For magnetized matter, the total force on an element of volume V is 
then 

ftot - / jtot X BdF + / (M X n) X Bd^ (46) 

Jv JdV 

where (neglecting the displacement current) jtot = j + V x M. Inserting 
the expression for the total current into the volume integral and using the 
divergence theorem on the surface integral, we obtain the force density 

Ftot - j X B + MfeVS^ (47) 

identical to the lowest order description from the Pauli equation (see Eq. 
(43)). Inserting the free current expression (45), due to Ampere's law, we 
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can write the total force density according to 

F^ ^-d^(—-M-B] +dk(H'B''). (48) 

V2mo / 

The first gradient term in Eq. (48) can be interpreted as the force due to 
a potential (the energy of the magnetic field and the magnetization vector 
in that field) , while the second divergence term is the anisotropic magnetic 
pressure effect. Noting that the spatial part of the stress tensor takes the 
form^^ 

T''' = -H'B'' + {B^/2na-M-B)S''', (49) 

we see that the total force density on the magnetized fluid element can be 
written F* = —dkT^'^, as expected. Thus, the Pauli theory results in the 
same type of conservation laws as the macroscopic theory. The momentum 
conservation equation (43) then reads 

p(|+V.v)v^-v(^-M.B)+5'^a.H-Vp, (50) 

where for the sake of clarity we have assumed an isotropic pressure, dropped 
the displacement current term in accordance with the nonrelativistic as- 
sumption, and neglected the Bohm potential (these terms can of course 
simply be added to (50)). This concludes the discussion of the spin-MHD 
plasma case. Next, we will look at some applications of the derived equa- 
tions. However, it should be noted that in many cases the spins are close 
to thermodynamic equilibrium, and we can thus write the paramagnetic 
electron response in terms of the magnetization^^ 

M = ^ tanh (^ B, (51) 

instead of using the full spin dynamics. Here B denotes the magnitude of 
the magnetic field and B is a unit vector in the direction of the magnetic 
field, ks is Boltzmann's constant, and T is the electron temperature. 



5. Examples and applications 

The above equations are quite complicated, but as such also extremely rich. 
Suitable and physically relevant approximations, such as the magnetization 
given by the expression (51), will however lead to considerable simplifi- 
cations. Below we consider two specific examples where such simplifying 
assumptions lead to interesting spin effects. 
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5.1. Spin solitons 

In Ref. 21, it was shown that the electron spin can introduce novel nonlinear 
structures in plasmas, with no limiting classical counterpart. In particular, 
the MHD limit for a electron-positron pair plasma is considered. 

Neglecting dissipative effect, the governing equations for the system of 
interest read 



dp 
at 



V • (pv) = 0, 



(52) 



/o(^+V-V)V=jxB-V-n 






2to2 



■ spin? 



(53) 



where fl = [(Tg + Tp)/2me]\ + {m? / p)'} (81 j is the total pressure tensor in 
the centre-of-mass frame and 



atanhf^^VB, 



(54) 



where we have assumed equal temperature T of the electrons and positrons. 
Moreover, we have 



dB 

'dt 

while the current is given by 



V X (V X B) , 



(55) 



— VxB-eo^(VxB)-^Vx 
po ot 2m 



.t.„h(gf)B- 



(56) 



For one-dimensional Alfvcn waves, the above system can be reduces to the 
modified Korteweg-de Vries equation^ ^ 



d 

^+CAsp 



c\ cos^ 9 d^ 






Vy 



[PBB^f v^ dVy 

2mekBT c\ d^ 



(57) 



where ca = [c^i?o/(c^MoPo + -Bq)] is the Alfvcn speed, c^.sp = c\/{l + 
Ssp) is the spin-modified Alfven speed, lOc = eBa/me is the cyclotron fre- 
quency, Bq is the magnitude of the unperturbed magnetic field, po is the 
unperturbed density, and we have the spatial coordinate ^ — x sin 0+z cos 0. 
We see that neglecting the spin contribution leads to a purely dispersive 
equations. Thus, the spin enables the formation of solitons with no limiting 
classical solution. 
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5.2. Ferromagnetic plasma behaviour 

For an ion-clcctron plasma, wc have the governing equations^^ 



^ + V-(pV)^0, (58) 



the momentum equation 



i + y.v)y = -v(|i-M.B)-Vp 

VMo / 2mem., \ y/p J 

and the ideahzed Ohm's law 

— ^VxiVxB), (60) 

where the variables are defined as above. Using the magnetization (51), we 
obtain a closed set of equations. 

In what follows, we will study the linear modes of this system, with a 
particular focus on the stability properties. With p = po+pi, B = Bq + Bi, 
M = Mo + Ml, and v = Vi, such that pi < po, l-Bil < l-Bol, 
|Mi| ^ |Mo|, and B^ = BqZ, we linearize our equations in the perturbed 
variables. Assuming that the background quantities are constants, the gen- 
eral dispersion relation can, after a Fourier decomposition, be written 



{-' ~ klcl) 



u^ - eel - klVlik)) (c.2 - klVlik)) + klklvXik) 



where Ca is the spin-modified Alfven velocity given by 





^A — 


1 + {nwlJ2m,c^J^})t&nh{iiBBQ/kBT) 


1/2 


Ca is 


the standard Alfven velocity Ca ~ {Bq/hoPoY^^ , 




Viik) = Viik) tanh -——- , 

TTli \ ksl J 


and 













= 
(61) 

(62) 



(63) 



(64) 



Here tOpf, — (pqc'^ /eomeirii)^''^ is the plasma frequency, uice is the electron 
cyclotron frequency associated with the external magnetic field (i.e. with the 
contribution to Bq from the spin sources excluded). The relation between 
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the full electron cyclotron frequency LOce ~ eBf)/ra^ and ojce is given by 
uJce = ^ce + {hu!pg/mec'^)teinh{fj,BBo/kBT). We stress that Va, which to 
some extent can be considered as an effective acoustic velocity, may be 
imaginary for a strongly magnetized plasma due to the spin contribution, 
a fact which will be explored in some detail below. 

We consider propagation perpendicular to the external magnetic field, 
which is the geometry which leads to instability most easily. For the case 
k — kx, Eq. (61) reduces to 



UJ 



(-12 



1 + {hu;^j2meC^uji"J)ta.nh{nBBo/kBT) 
Amirrie mi \ kBT 



(65) 



The necessary and sufficient instability condition can thus be written as 
Psp + Pm + P + y^ < 0, (66) 

where the total pressure Ptot = Psp + Pm + P consists of the effective spin 
pressure Pgp = —{pahajce/ini) ta,nh{pB Bq / kBT) , which is the only negative 
pressure term and therefore the source of the instability, the magnetic pres- 
sure Pin = poC\/[l + {hL0p^/2mec'^i0ce )tSi''nh{fiBBo/kBT)], and the parti- 
cle pressure P = noniicl, containing both the thermal and Fermi pressure 
part. Thus, a plasma can contain a magnetization instability, much like the 
gravitationally induced Jeans instability.^^ 



6. Conclusions and future possibilities 

We have seen that quantum effects, and in particular the electron spin, can 
introduce new and interesting aspects in plasma theory and experiments. 
It is expected that the rapid development of quantum plasma theory will 
be fueled by recent experiments (see, e.g, Ref. 31), and that new regimes 
of interest will enter the arena as the next-generation laser systems gets 
online. 

Some particular developments that would be of interest is to further 
look at spin effects from a kinetic perspective. ®^'^^ Such treatments would 
be similar to the density matrix approach, and could, by using analogies 
from classical kinetic theory, spur the experimental interest in collective 
quantum effects and the transition from quantum to classical behaviour. 
Moreover, the fluid equations presented here has not been analyzed to their 
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full extent. For example, there are terms which have been neglected, that 
could produce interesting nonlinear effects in quantum plasmas, such as 
spin self-interaction. Such topics will be approached in future research. 
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